Magnetic properties of pure and Gd doped EuO probed by NMR 
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An Eu NMR study in the ferromagnetic phase of pure and Gd doped EuO was performed. A 
complete description of the NMR lineshape of pure EuO allowed for the influence of doping EuO 
with Gd impurities to be highlighted. The presence of a temperature dependent static magnetic 
inhomogeneity in Gd doped EuO was demonstrated by studying the temperature dependence of 
the lineshapes. The results suggest that the inhomogeneity in 0.6% Gd doped EuO is linked to 
colossal magnetoresistance. The measurement of the spin-lattice relaxation times as a function of 
temperature led to the determination of the value of the exchange integral J as a function of Gd 
doping. It was found that J is temperature independent and spatially homogeneous for all the 
samples and that its value increases abruptly with increasing Gd doping. 

PACS numbers: 71.30.+h, 75.30.m, 76.60.k, 75.50.Pp 



I. INTRODUCTION 

Numerous studies on EuO have been published since 
the discovery of the ferromagnetic phase of this com- 
pound by Matthias et al. in 1970," and still very re- 
cently several experimental research projects have been 
conducted on pure and Gd doped EuO j 2 ' 3 i 4 i 5 There are 
several reasons for this interest: first, EuO is one of the 
few natural ferromagnetic semiconductors. Second, there 
is currently a great deal of attention on ferromagnetic 
semiconductors like (Ga,Mn)As and (In,Mn)As. 6 Third, 
a new field called spintronics has been developed around 
the possibilities of using the spin degree of freedom of 
the electron in solid-state electronics Fourth, EuO is 
an ideal system for testing new theories in magnetism, in 
particular the recent developments made on the Kondo- 
lattice model. The localized magnetic moments of the 
half-filled 4/-shell of the Eu atoms and the existence of 
a conduction band makes EuO an appropriate system to 
test the Kondo-lattice modeli&Siifli In addition, the low 
magnetic anisotropy of the material along with the local- 
ized J — S = 7/2 spins of the Eu 2+ ions makes europium 
monoxide a nearly ideal Heisenberg ferromagnet. 

Another peculiarity of EuO compounds such as O- 
rich EuO and Gd doped EuO is their colossal mag- 
netoresistance (CMR)i"i 12 i 13 i 14 i 15 CMR has also been 
observed in manganites. Nuclear magnetic resonance 
(NMR) has been proven to be a good technique to study 
local magnetic microstructures such as stripes and mag- 
netic polarons>iSiii However, to our knowledge, such a 



study in europium chalcogenides has never been per- 
formed and NMR data in EuO have been published only 
for temperatures far below the transition temperature, 
i.e. far below the temperatures at which CMR effect is 
observed in electron doped EuO. The main reason that 
led us to study europium chalcogenides rather than man- 
ganites is that the former are simple diatomic cubic crys- 
tals whereas the structure of the latter is substantially 
more complex. Consequently, studies of the (Eu,Gd)0 
system should allow us to extract detailed information 
from NMR measurements without the complications of 
numerous exotic phenomena of manganites such as phase 
separation, charge and orbital ordering, and Jahn- Teller 
distortions. 



In this paper we present NMR results on single crystals 
of pure and Gd doped EuO. In Gd doped EuO, Gd atoms 
play the role of electron donors and are expected to affect 
the localized magnetism of EuO negligibly. One goal of 
this study was to examine the effect of Gd doping on the 
magnetic properties of EuO. In SectionlLTI we describe 
the characteristics of the samples we studied and give a 
detailed analysis of the NMR lineshape of Eu in EuO. In 
Seclion lTTTI we describe the temperature dependence of 
the relaxation times of the Eu nuclear spins. Finally, in 
Section lTVI we discuss the properties of the lineshape of 
Eu in Gd doped EuO and discuss some models to link 
our observations to CMR behavior. 
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II. SAMPLE CHARACTERIZATION 

We studied four types of samples: pure EuO, 0.6%, 2%, 
and 4.3% Gd doped EuO. According to Samokhvalov et 
al, these samples can be categorized as follows: 14 samples 
of pure and stoichiometric EuO arc insulators; samples 
with a Gd concentration smaller than about 1.5% un- 
dergo a metal-insulator transition (MIT) when the tem- 
perature is increased above about 30 K, the low temper- 
ature regime being metallic; samples containing a Gd 
concentration larger than about 1.5% have a metallic be- 
havior at all temperatures. Godart et al. also observe 
similar behavior in Gd doped EuOj 1 ^ but not Schoenes 
and Wachtcr^ 



A. Spectrometer 

The spin-spin relaxation time T2 of 151 Eu and 153 Eu 
nuclei in pure and Gd doped EuO can be very short, es- 
pecially when the temperature is increased towards the 
transition temperature. To measure signals with short 
T2, we built a spectrometer with the following specifics: 
minimum pulse length of 20 ns, minimum delay between 
two pulses of 300 ns, recovery time of the receiver of 
300 ns, and maximum sampling rate of the receiver of 
lGSample/s. The high sampling rate is needed because 
of the short echoes^ 

In order to detect the signal using a short delay be- 
tween the two excitation pulses of a spin-echo sequence, 
it was necessary to use a low Q tank circuit. Typically, 
we used a Q of about 10, which corresponds to a recov- 
ery time of about 300 ns. Note that the use of such a 
small value of Q is possible because of the presence of 
an amplification factor, a property inherent to magnetic 
materials. 



B. Remarks on amplification factor 

In ferromagnetic materials, the excitation RF field Hi 
acting on the nuclear spins is amplified by a factor 77 
via the magnetic susceptibility of the unpaired electron 
spinsiSS In EuO the magnetic anisotropy is low and con- 
sequently the amplification factor for nuclei in domain 
walls is not substantially stronger than for nuclei in do- 
mains. Indeed, the amplification in domains is given by 

^Domain — 

\H hf \/\H an + H mt \, where H hf is the hy- 
perfine field, H an is the anisotropy field and Hi nt is the 
internal field defined as the sum of the external field, 
the demagnetization field and the Lorentz field 2 *!-. The 
amplification in a domain wall of a spherical particle is 
given by r] DW = nD\H hf \/(N5\M(T)\), where N is the 
demagnetization factor of the domain, D is the domain 
size, 6 is the width of the domain wall and M(T) is the 
magnetization in domains at temperature Tj2£ and there- 



fore: 

VDomain _N 5 \M{T)\ 

r\ DW 7T D \H an \ 

From the study of the dynamic susceptibility on an 
EuO sphere, Flosdorff et al. have determined that 
S/D 0.08|M(T = 0)\/\M(T)\M Then, by taking 
N = 47r/3, assuming a spherical domain, \M(T 
1.2 K = 1623 emu/cm 3 ; 24 and \H an (T = 4.2 K)| = 
247.5 ± 2.5 Oe^ we obtain r]Domain/r]DW — 0.7. There- 
fore, at 4.2 K the amplification in domains is only slightly 
weaker than the amplification in domain walls. More- 
over, since \H an (T)\ cx (\M(T)\/\M(T = 0)|) 5 ^ the 
ratio increases rather rapidly with increasing temper- 
ature. It appears then that we cannot a priori decide if 
the NMR signal will be dominated by nuclei in domains 
or nuclei in domain walls. However, several considera- 
tions indicate that what we observe does not depend on 
whether the nuclei are located in walls or in domain walls. 
First, the lineshape we observed is similar to the one ob- 
served in previous experiments done in field high enough 
to suppress domain walls i2L22i22i Second, in a field of 4 T, 
when all domain walls are suppressed, we observed the 
same temperature dependence of the relaxation times as 
in zero-field. Therefore, in zero-field we did not observe 
any additional relaxation mechanism related to the pres- 
ence of domain walls such as a relaxation mechanism due 
to domain walls motion. 



C. Magnetization vs. temperature 

In ferromagnetic materials, the hyperfine field is pro- 
portional to the magnetization. Therefore, if we assume 
that the nuclear dipolar field is negligible compared to 
the hyperfine field, a measurement of the temperature 
dependence of the zero-field 153 Eu NMR frequency gives 
the temperature dependence of the magnetization. The 
frequency measurements are presented in Fig.^ We de- 
termined the Curie temperature by measuring the sus- 
ceptibility x m the paramagnetic regime, plotting 1/x as 
a function of temperature and determining the temper- 
ature at which the extrapolated 1/x line reaches zero. 
This temperature corresponds to the so-called paramag- 
netic Curie temperature 9c- It is important to note that, 
according to the results of Samokhvalov^ in Gd doped 
EuO samples the value of 6*c is larger than the value of 
the Curie temperature Tc, the temperature above which 
the spontaneous magnetization vanishes. In particular, 
the value of Tc of both pure and 0.6% Gd doped EuO is 
about 69.55 K. 

Although the magnetization in these systems is not 
well described by mean-field theory since spin waves 
play a crucial role at low temperatures, we plotted Bril- 
louin curves along with the data as rough approximations 
of the entire magnetization curves intended to serve as 
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Pure EuO (e c =69.55K) 
0.6% Gd doped EuO (9 C =73.05K) 
2% Gd doped EuO (e c =123.45K) 
4.3% Gd doped EuO (9 C =132.55K) 
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138 139 
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FIG. 1: Zero-field NMR frequency of the center of the NMR 
line of pure and Gd doped EuO vs. temperature. Brillouin 
functions are plotted along with the data. 

guides for the eyes. The data shown in Fig.Q] are com- 
parable to the data published by Mauger et alM> This 
is one of the indications that our samples are consistent 
with those of other authors, though coming from differ- 
ent sources at different times. Other such signs of con- 
sistency among samples will come from NMR data, as 
shown below. 

Kapusta et al. measured the temperature dependence 
of the frequency in various manganites and observed that 
the hyperfine field did not vanish at T = The au- 

thors concluded that there was a residual magnetization 
above Xc that could be due to the presence of ferromag- 
netic regions such as magnetic polarons. We could not 
observe the resonance near Xc, but our data differ qual- 
itatively since these authors observe nearly no drop in 
magnetization up to Tq whereas we observed that the 
resonance frequencies already decrease substantially if T 
is increased up to kTc- 

D. Lineshape in pure EuO at 4.2 K 

The NMR lineshape of 153 Eu in single crystals of pure 
EuO at 4.2 K has already been studied several times. In 
1966, Boyd found a single sharp resonance using con- 
tinuous wave NMRi^i Later, Raj et al. observed that 
the zero-field lineshape is composed of a sharp central 
line and two wings but they did not discuss the broad- 
ening mechanisms^ Then, Fekete et al. observed five 
quadrupolar lines in a saturation field of 2T^£ Finally, 
Arons et al. observed a zero-field lineshape that they de- 
scribed as a sharp line on top of a broad linem&22& They 
claimed that the sharp line corresponded to the signal 
from nuclei located in domain walls and that the broad 
line corresponded to the signal from nuclei in domains. 



FIG. 2: Zero-field lineshape of 15d Eu in EuO at 4.2 K: the 
black dots correspond to the FT of the echo and the open 
circles to the point-by-point measurement. The dotted line is 
the fit of the data using © and the plain line is the same fit 
corrected for Tb effects. 



We present in Fig.[2]the zero-field lineshape of 153 Eu 
in pure EuO that we measured at 4.2 K. We determined 
the lineshape by two different techniques: on one hand, 
we measured the echo integral at a certain number of dis- 
crete values of frequency (point-by-point measurement). 
Alternatively, we performed a Fourier transform (FT) of 
the echo measured at the frequency of the central line. 
As shown in Fig.[3 both methods give an almost identical 
lineshape. 

Since contradictory results and analysis of the line- 
shape of 153 Eu in pure EuO at 4.2 K were published, 
we present here our own analysis. Because the Eu sites 
in EuO have a nominal cubic symmetry, we expected 
to observe a single narrow line. However, as shown in 
Fig.[3 we observed an intense central peak and two wings 
with a structure suggesting that each wing is composed 
of two broadened lines. In order to distinguish between 
the broadening due to electric field gradients (EFGs) and 
the broadening due to the magnetic environment of the 
nuclei, we took advantage of the fact that there are two 
isotopes of europium with similar natural abundance. 
The gyromagnetic ratio and the electric quadrupole mo- 
ment of the two isotopes are shown in TableQ] The zero- 
field lineshape of 151 Eu in pure EuO at 4.2 K is shown 
in Fig. [31 By comparing the 153 Eu and the 151 Eu line 
at 4.2 K we observed that the lines cannot be described 
solely by magnetic broadening or by quadrupole electrical 
broadening. We therefore came to the conclusion that the 
broadening was due to a combination of both effects. We 
will discuss the origin of such a broadening in Sect. II V Al 

In the presence of EFGs, since 153 Eu has a spin / = 
5/2, we expected to observe five lines all separated by the 
same frequency interval Avq and with intensity ratios 
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TABLE I: Spin, natural abundance, gyromagnetic ratio, and 
electric quadrupole moment of the two europium isotopes^ 

Spin Nat. abund. (%) 7„/27r (MHz/T) Q (barn) 



^EuT 

151 Eu 



5/2 
5/2 



52.19 
47.81 



4.6745 
10.5856 



2.41 
0.903 



E 
< 




310 



312 

Frequency (MHz) 



314 



5:8:9:8:5i2£ We tried to reproduce the observed lineshape 
theoretically by assuming a distribution in EFGs and a 
magnetic broadening. We described the inhomogeneities 
in the EFGs with a Gaussian distribution leading to a 
distribution in frequency G{v) of second moment Sau q 
that broadens the five lines. As will be discussed in Sec- 



tion llll Bl a Lorcntzian distribution describes correctly 
the magnetic broadening. The broadening resulting from 
these two mechanisms is a convolution of the two distri- 
bution functions. The convolution of a Gaussian and a 
Lorentzian gives a Voigt function, i.e., a function of the 
form 



V{v) 



7r 3 / 2 



(v - v Q - V2at) 2 + ary 



■dt, (2) 



FIG. 3: Zero-field lineshape of 151 Eu in EuO at 4.2 K: the 
black dots correspond to the point-by-point measurement. 
The dotted line is the computed lineshape and the plain line 
is the same computed lineshape corrected for T2 effects. 



J 



where T is the full width at half maximum (FWHM) 
of the Lorentzian and a is the second moment of the 
Gaussian^ The resulting lineshape is then given by 



/(") = 



m=2 9 1 

m=-2 



iy - (>o 



where T represents the FWHM of the magnetic distribu- 
tion and the factor (1 — m 2 /9) accounts for the differ- 
ence in intensity of the five quadrupolar lines. By fitting 
the measured 153 Eu lineshape with this function, it was 
possible to determine the parameters T, Aj/q and Sav q ■ 
The fit is represented by a dotted line in Fig.[3 and the 
deduced values of the fitting parameters are shown in 
Table ITU However, there was another characteristic that 
needed to be taken into account: Raj et al. showed that 
the Eu spin-spin relaxation time in pure EuO is short 
and frequency dependent and, as a consequence, some 
of the nuclei might not be observed if the delay between 
the excitation pulses is too long. 3 - 4 . They observed that 
the intensity of the central part of the line is strongly re- 
duced. We demonstrated this phenomena in Fig.0]where 
we plotted the FT of the echo measured with two different 
delays*! 3 - We took this effect into account in the compu- 
tation of the lineshape of 153 Eu by multiplying 10 by an 
approximated shape of the frequency distribution of T2 





-t 2 
e 1 


(3) 


- mAi/g) - 


- V2\m\6 AvQ tf + (|r)2 ' 


1 

TABLE II 
153 Eu alon; 


Values of the lineshape fitting 
5 with the deduced parameters for 


parameters for 
151 Eu. 




T [MHz] Auq [MHz] 


8 Avo [MHz] 


lby Eu 
151 Eu 


0.072 0.75 
0.163 0.281 


0.2 
0.075 



at 4.2 K deduced from our measurements. We will discuss 
the physical reason for this T2 distribution in Sect. lIH 51 
The computed lineshape is shown in plain line in Fig.[21 
along with the measured lineshape. 

To verify if our description of the 153 Eu lineshape was 
correct, we computed the 151 Eu lineshape from the fitting 
parameters that we determined for 153 Eu and compared 
the resulting curve with the measured 151 Eu lineshape. 
We first calculated the parameters for 151 Eu using the 
relations 
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FIG. 4: Zero-field lineshape of 153 Eu in EuO at 4.2 K ob- 
tained by FT. The black dots correspond to a delay of 10.5 fis 
and the open circles to a delay of 15.5 /xs. This shows the pres- 
ence of frequency dependent T2's (the Tb's are shorter near 
the line center). 



Ti51 



In 



153 
151/ 



In 



ri53, 



153 Q 
151 Q 



153, 



153?5 °Af<3,153, 



(4) 
(5) 
(6) 



where the indices 151 and 153 refer to 151 Eu and 153 Eu 
respectively. The computed lineshape is shown as a 
dotted line in Fig.EI Then, we deduced the T 2 distri- 
bution of 151 Eu from the T2 distribution of 153 Eu us- 
ing the following process: first, since we observed that 
T2 at 4.2 K was inversely proportional to 7^ (c.f. Sec- 
tio nllll B]l . we divided the amplitude of the distribution 
by ( 151 7„/ 153 7„) 2 . Second, since T2 depends on j n (and 
is therefore of magnetic origin), we expected the T2 dis- 
tribution to scale with the magnetic distribution. We 
thus calculated the 151 Eu T 2 distribution by multiplying 
the width of the 153 Eu T 2 distribution by 151 7„/ 153 7„. 
The deduced distribution is in good agreement with our 
data. In order to take into account the fact that the 
T2 of the central peak was too short to be detected, we 
introduced a cutoff in the computed distribution curve. 
Finally, the T 2 corrected 151 Eu lineshape was obtained by 
multiplying the computed lineshape by the resulting T 2 
distribution. The result is plotted in solid line in Fig.|3| 
along with the measured data. The calculated line was 
in reasonable agreement with the data, thus confirming 
our description of the lineshape in terms of a magnetic 
and a quadrupolar broadening^ 

The origin of the quadrupolar splitting is most likely 



intrinsic to EuO and not due to oxygen vacancies as sug- 
gested by Arons et al£& Oxygen vacancies will lead to a 
clear breaking of the cubic symmetry and therefore, since 
the quadrupolar moment of Eu is large, the quadrupo- 
lar splitting is expected to be large and the Eu sites lo- 
cated next to the vacancies will not contribute to the ob- 
served line. Note that the electric quadrupole splitting 
we observed in zero-field is similar to the one observed by 
Fekete et al. in an external field Hq along the easy axis 
[111] saturating the magnetization (2T< \H \ <3T)i21 



III. EFFECTS OF MAGNONS AND DOPING 

The measurement of NMR relaxation times in eu- 
ropium chalcogenides gives some valuable information 
about the dynamics of their electron spins through the 
hyperfine coupling. Boyd was the first to report NMR 
data on EuS^ and together with Charap, used spin- 
wave theory to deduce values of Jj and J 2 , respec- 
tively the nearest neighbor and next nearest neighbor 
exchange constant, from the temperature behavior of the 
NMR frequency of both 151 Eu and 153 Euii The first val- 
ues of NMR relaxation times in EuO were published in 
1965 by Uriano and Streever. 42 Following this work, sev- 
eral groups studied NMR relaxation times and compared 
their data to spin- wave theory and to the Suhl-Nakamura 
theory in order to explain their results and to obtain in- 
formation on the magnetic properties of the system. All 
the results published concern measurements performed at 
T < 20.3 K. We extended the study of relaxation times 
to much higher temperatures and to Gd doped EuO. We 
determined the influence of doping on the dynamic prop- 
erties of the electronic magnetization by comparing the 
relaxation times in EuO to the relaxation times in Gd 
doped EuO. 



A. Spin-lattice relaxation times 

We used a saturation recovery sequence to perform Ti 
measurements. The sequence destroys the nuclear mag- 
netization by applying a series of pulses before measuring 
the signal. It consists of a series of pulses of 0.2/is sepa- 
rated by a time of the order of T 2 forcing the nuclear spins 
to loose their coherence. A standard spin-echo sequence 
is applied after a while to detect the amplitude of the 
magnetization that has recovered equilibrium. Note that, 
unlike for T 2 measurements, no substantial frequency de- 
pendence of Ti was observed, neither in pure EuO nor in 
Gd doped EuO. 

Barak, Gabai and Kaplan determined that the spin- 
lattice relaxation in powdered EuO in an external mag- 
netic field at temperatures below 14 K are dominated by 
two-magnon processes^ 3 - Also, between 14 K and 20 K, 
they observed that a three-magnon process surpasses the 
two-magnon process and that it is enhanced by a second- 
order three-magnon process. The results we obtained at 
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FIG. 5: Zero-field spin-lattice relaxation rates of Eu in 
pure and Gd doped EuO as a function of temperature. The 
expression used to fit the curves is defined in JjJ. 



zero-field confirmed that, for temperatures above about 
14 K, the three-magnon process is the dominant spin- 
lattice relaxation process in EuO. The relaxation rate 
(Fig.|SJ was found proportional to T 7 / 2 . This is con- 
sistent with the temperature dependence of the theoret- 
ical expression first derived by Oguchi and Keffer for a 
three-magnon process44 We rederived the expression and 
adapted it to the case of zero-field measurements to ob- 
tain (Appendix 0: 

1 _ 11-29 A 2 fk B T\ 7/2 
Ti 16(2tt) 5 2 JS ■ HS \ 2JS J ' U 

where A is the hyperfine constant, J = Ji + J 2 the ex- 
change integral, S the electron spin, and k B the Boltz- 
mann constant. As for the results of Barak et at, we 
had to multiply the expression by the exchange scat- 
tering enhancement factor that we determined to be 
£ = 8, in agreement with the calculation of Bccman and 
Pincus. 45 Taking S = 7/2 and determining A from the 
zero-field resonance frequency extrapolated to T = OK, 
A = -hw/S = -2.6 • 10~ 26 J, we deduced from the fit a 
value of the exchange integral of J/k B — 0.755 ± 0.01 K. 

NMR thus provides an accurate determination of J. 
Even if the enhancement factor is not precisely known, 
the fact that 1/Ti, and therefore £/Ti, is proportional to 
j-9/2 ^ c £ 0) means that a large change in £ leads to 
a small change in J. The influence of a variation in the 
enhancement factor on the determination of J is there- 
fore limited, and we evaluated the error to be less than 
1%. The other causes of uncertainty are the error on 
the Ti measurement and the error on the temperature 
measurement. An evaluation of the total error leads to 
AJ/k B = 0.01 K. 



TABLE III: Values of J as a function of Gd doping. 



Doping level x (%) 


Exchange integral J/fcs (K) 





0.755 ±0.01 


0.6 


0.750 ± 0.01 


2 


1.205 ±0.01 


4.3 


1.210 ±0.01 



Previous measurements on powdered EuO by neutron 
scattering experiments and specific heat measurements 
led respectively to (Ji + J 2 )/k B = 0.725 ± 0.006 
and (Ji + J 2 )/k B = 0.714 ± 0.007K. 47 More recently, 
Mook measured single crystals of EuO by neutron scat- 
tering methods and he obtained the following values: 
J x /k B = 0.625 ±0.007K and J 2 /k B = 0.125 ± 0.01 K. 48 
These values lead to a value of the exchange integral 
J = 0.750 ± 0.017 K which is in very good agreement 
with our measurement. Note that in 1966 E. L. Boyd 
deduced J\ and J 2 by measuring the temperature de- 
pendence of the NMR frequency of 153 Eu>2£ He found a 
value of Ji/k B = 0.750 ± 0.0025 K that is very close to 
the value of J that we determined. However, he deduced 
a negative value of J 2 . 

In Gd doped EuO (Fig.EJl, we observed that the tem- 
perature dependence of the spin-lattice relaxation times 
is the same T 7 / 2 . Hence, as for pure EuO, the three- 
magnon process is the dominant relaxation process. Us- 
ing expression Q and the enhancement factor £ — 8, we 
deduced the values of J shown in Table iHTl for the doped 
samples. Clearly, the samples can be separated in two 
categories according to their value of J (Fig.[5J| . 

To the best of our knowledge, no direct measurement 
of the exchange integral in Gd doped EuO has ever been 
published and our results provide therefore the first mea- 
surement of J in Gd doped EuO. 



B. Spin-spin relaxation times 

We performed spin-spin relaxation measurements in 
pure EuO and in Gd doped EuO using a standard spin- 
echo sequence composed of two consecutive pulses of du- 
ration t\ and t 2 separated by a varying delay. Typical 
values were t\ = 0.1 /xs and t 2 = 0.2 /is. Since T 2 is 
strongly frequency dependent, it is necessary to specify 
at what frequency it was measured. In the following, 
T 2 always corresponds to the spin-spin relaxation time 
measured at the frequency corresponding to the center 
of the central line, or the center of gravity of the line in 
case the central line is not well defined. Note that all 
the relaxation times that will be presented in the follow- 
ing correspond to measurements done with short delays. 
Therefore, the discussion and the analysis will focus on 
the fast decay times unless specified otherwise. 

Except for the value of T 2 at 20.3 K determined by 
Uriano and Streever— and the decay curve at 13.8 K 
published by Barak et al.^ there are no spin-spin re- 
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laxation time data for EuO above 4.2 K in the literature. 
However, several studies of T 2 at 4.2 K and lower tem- 
peratures were performed following the work of Barak 
et aZ.i 49 i 50 Raj et aZ.jSi Fekete et ai.^i the Liitgemeier 
group; 28 i 29 i 39 i 52 i 53 and more recently Pieper et alJ& The 
main reason for this interest is that the decay curve is not 
a single exponential. As Liitgemeier et alM first noted, 
this is because there are two relaxation mechanisms: a 
fast one due to the Suhl-Nakamura indirect interaction, 
and a slow one due to direct dipolar coupling. However, 
the original Suhl-Nakamura theory fails to explain the 
observed spin-spin relaxation time and its frequency de- 
pendence that was first reported by Raj et al.^ Barak et 
al. used the theory developed by Hone, Jaccarino, Ngwe 
and Pincus^(in the following we will refer to this theory 
as the HJNP theory) to explain why the value of T 2 they 
observed in a powdered EuO sample at 4.2 K is longer 
than the one predicted by the Suhl-Nakamura theory^ 9 . 
The HJNP theory assumes the existence of an inhomoge- 
neous line broadening resulting from random microscopic 
inhomogeneities and predicts a frequency dependent T 2 . 
The HJNP theory predicts a Lorentzian line shape and 
for this reason we used a Lorentzian as magnetic distribu- 
tion function in the calculation of the lineshape presented 
in Section llTDl 

In 1976, Fekete et al. successfully described the de- 
cay curves of the five transitions they had observed by 
assuming an inhomogeneous Suhl-Nakamura relaxation 
process (deduced from the HPJN model) and a dipolar 
coupling.— They took into account all the elements nec- 
essary to describe the lineshape and its variation with 
delay as shown in Fig.0J However, they did not discuss 
the implication of their results on the lineshape. Arons 
et al. published in 1975 the zero-field lineshape vs. delay 
and observed, as we did, the reduction of the amplitude of 
the central peak with increasing delayi 2 ^ However, Arons 
et al. claimed that the T 2 of the central line was shorter 
because the signal was coming from nuclei located in do- 
main walls. Our combined results on T\ and T 2 mea- 
surements as a function of temperature point to more 
intrinsic mechanisms. 

We measured the temperature dependence of the spin- 
spin relaxation times (Fig.|B|). For pure EuO as well as 
for all the Gd doped samples, T 2 seems to be tempera- 
ture independent at temperatures below about 15 K. For 
higher temperatures, we observed a rapid increase of the 
relaxation rates with increasing temperature. We ana- 
lyze these two regimes below. 



nuclear spins, in particular by processes involving three 
magnons. One way to describe these processes is to con- 
sider t he fluctuations of the hyperfinc field. We can writ e 
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FIG. 6: Zero-field spin-spin relaxation rates of 153 Eu in pure 



and Gd doped EuO as a function of temperature, 
rameter a and the fitting curve are defined in @ . 



the perturbation Hamiltonian as 



The pa- 



Hlf{t) = 



- ln h[H x {t)i x + Hy(t)iy + H z (t)l Z 



(8) 



where H x (t), H y (t) and H z (t) are the fluctuating fields 
at the nuclear site due to fluctuations of S in the x, 
y and z directions respectively^ The relaxation rates 
deriving from these interactions can be calculated using 
the Redfield theory^ The calculation for the particular 
case of a spin / = | gives the following result for the 
spin-spin relaxation time: 



1 2T72 , 1 

J 2 Jl 



(9) 



1. Lifetime effect 

We saw in Scct. lIII Al that the spin-lattice relaxation 
processes are dominated by the scattering of magnons by 



where H | is the amplitude of the correlation function be- 
tween H z (t) and H z (t + At), To is the correlation time (or 
lifetime) of the scattering process between the magnons 
and the nuclear spins, and 
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5 for the —5/2 «-> —3/2 and 3/2 «-> 5/2 transitions 
= 7(J + 1) - m(m + 1) = <( 8 for the -3/2 <-> -1/2 and 1/2 <-> 3/2 transitions (10) 

9 for the — 1/2 <-» 1/2 transition. 



We have assumed that the fluctuations of the three com- 
ponents of field x, y, and z are independent and that 
the correlation functions are simple exponential. 37 Since 
the length of the pulses as well as the amplitude of the 
RF field were chosen such as to excite all the transitions, 
we expect to have a mixture of these rates. As shown 
in Fig.El t ne factor a obtained from fitting the curves 
with the function I/T2 = a/T\ + (3 is between 6.3 and 
9.9. This is in good agreement with the theoretical val- 
ues (|1(J[1. Therefore, according to the Redfield theory, 
the temperature dependence of the spin-spin relaxation 
in pure and Gd doped EuO is entirely determined by 
transverse fluctuating fields. Our results showed that 
these fluctuations were due to fluctuations of the elec- 
tronic spins well described by spin-wave theory. 



2. Temperature independent mechanisms 

For temperatures far below the magnetic transition 
temperature, the Suhl-Nakamura interaction is expected 
to be temperature independent since its temperature de- 
pendence comes mostly from the hyperfine constant 
The nuclear dipole-dipole interaction is also independent 
of temperature. Therefore, two processes compete and 
may be the source of the observed temperature inde- 
pendent T2S. In the case of pure EuO, the relaxation 
time due to the dipole-dipole interaction is considerably 
slower and we observed, in agreement with previous mea- 
surements, that the Suhl-Nukamura processes dominate 
the dipole-dipole processes for short delays. As shown in 
Fig. El the relaxation time of 153 Eu in Gd doped EuO is 
substantially slower than in pure EuO. In 2% Gd doped 
EuO at 3.8 K, we have measured T2 = 183 fis. This value 
is very close to the value T2 = 194 ± 2.4 /j,s deduced, fol- 
lowing the treatment of Bohn et al. ^ from an exponen- 
tial fit of the slow spin-echo decay measured in pure EuO 
and thought to be due to dipole-dipole interactions*^ 

We observed that adding Gd in the EuO matrix leads 
to magnetic broadening, which increases with doping (c.f. 
Section EE). Therefore, in agreement with the HJNP 
theory, the replacement of Eu atoms by Gd atoms in 
the EuO matrix reduces the allowed mutual spin-flips 
between Eu nuclei leading to the Suhl-Nakamura relax- 
ation and consequently we expect the relaxation rates 
to be smaller when the concentration of Gd is increased. 
The relaxation through dipole-dipole interactions is how- 
ever still effective since in addition to the vanishing 1^17 

terms, the interaction Hamiltonian has diagonal I z ,ilg,j 
terms that do not vanish even if the frequencies of spins 



i and j are very different. Note that the temperature 
independent relaxation rate is larger for 4.3% Gd doped 
EuO than for 2% Gd doped EuO. The presence of an 
increasing number of conduction electrons when the Gd 
concentration is increased above about 1.5% might intro- 
duce an additional relaxation process through an RKKY 
indirect interaction. 



IV. MAGNETIC INHOMOGENEITIES 

A. Temperature dependence of the pure EuO 
lineshape 

To the best of our knowledge, no measurement of the 
lineshape of pure EuO above 4.2 K has ever been pub- 
lished. We report here the first study of the temperature 
dependence of the 153 Eu lineshape. In Fig.0 we present 
the lineshape of pure EuO at 4.2 K and at several tem- 
peratures between 4.2 K and 42 K. These lines were ob- 
tained by Fourier transforming the spin-echo measured at 
the frequency of the center of the lines. Since, as shown 
in Section lTl Gl the NMR frequency varies considerably 
with the temperature, we shifted the center of the lines 
presented in the figure to zero frequency in order to com- 
pare the lines recorded at different temperatures. We 
also adjusted the amplitude of the lines to superimpose 
the wings of all the lines. This adjustment allowed us to 
compare the amplitude and the width of the central peak 
of the lines. 

We observed that the intensity of the central peak was 
greatly reduced with increasing temperature. As noted in 
Section rill Bl the spin-spin relaxation times are frequency 
dependent and they are shorter for frequencies close to 
the central peak. Also, we observed that T% is strongly 
temperature dependent, decreasing with increasing tem- 
perature. Therefore, we concluded that the intensity of 
the central peak decreased with increasing temperature 
as a consequence of the decreasing value of T2, our spec- 
trometer not being able to detect all the nuclei of the 
central peak. We also observed that the width of the cen- 
tral peak increases almost linearly from about 0.11 MHz 
at 4.2 K to about 0.23 MHz at 41.4 K. This broadening 
mechanism did not seem to be linked to T2 effects since 
we measured T2 = 4/j.s at 41.4 K, which corresponds to 
a Lorentzian width of about 0.08 MHz. 

Finally, we observed that the electrical quadrupole 
splitting decreases with increasing temperature. We es- 
timated the change to be of the order of 25%, from 
Av Q S 0.75MHz at 4.2 K to Av Q = 0.55 MHz at 41.4K. 
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FIG. 7: Zero-field lineshape of 153 Eu in EuO vs. tempera- 
ture. The lineshapes were obtained by Fourier transform of 
the spin-echo signal. 
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FIG. 8: Zero-field lineshape of 153 Eu in 0.6% Gd doped EuO 
vs. temperature. The lineshapes were obtained by Fourier 
transform of the spin-echo signal. 



We think that this change might be due to a variation of 
the magnetostriction with temperature. This assumption 
is supported by the results of F. Levy who observed, by 
X-ray measurements, that the spontaneous magnetostric- 
tion of EuO decreases by about 25% between 4.2 K and 
40K. 58 Note that if the quadrupolar splitting was due to 
the presence of defects in the crystals, as was proposed 
by Arons et al.jEL we would probably not observe such a 
temperature dependence. 



B. Influence of Gd on EuO lineshape 

We now turn to the analysis of the lineshapes of 153 Eu 
in Gd doped EuO. We observed that the frequency of the 
central transition of the 0.6% Gd doped EuO lineshape at 
4.2 K was only slightly shifted towards higher frequency 
compared to pure EuO, from 138.45 MHz to 138.48MHz. 
However, the central line was substantially broader as 
can be seen on Fig.|S| As discussed in Section Ul Dl the 
width of the central transition is determined by mag- 
netic broadening only. We think that the observed line 
at 4.2 K is broader in the 0.6% sample than in the pure 
sample because of the random distribution of Gd atoms 
in the EuO matrix. Since different Eu atoms have differ- 
ent numbers of nearest and next nearest Gd neighbors, 
the local field acting on Eu atoms, which is influenced by 
the presence of Gd, is not the same at each Eu site. An 
other noticeable difference between pure EuO and 0.6% 
Gd doped EuO is that the quadrupolar structure is not 
resolvable in 0.6% Gd doped EuO. We think that this is 
also due to the random distribution of Gd atoms in EuO 
which leads to a random distribution of EFGs. 

We also measured the lineshape of 0.6% Gd doped EuO 
as a function of temperature from 4.2 K to 49.5 K (Fig.[SJ|. 



We observed a temperature dependent broadening of the 
central line. We adjusted the amplitude of the differ- 
ent lines by superimposing the wings of all the finest 
We plotted in Fig.HHI the full width at half maximum 
(FWHM) of the central line as a function of temperature 
for all the samples*^ It appears that for temperatures 
above about 30 K, the FWHM of the 0.6% Gd doped EuO 
lineshape increases sharply with increasing temperature. 
We did not observed this phenomena in pure EuO. We 
demonstrate that the broadening could not be explained 
by a variation in T2 and that the broadening was there- 
fore a static broadening. Indeed, the measured I/T2 is 
much less than the observed lincwidth. Moreover, we ob- 
served that T2 has the same temperature dependence in 
0.6% Gd doped EuO as in pure EuO (c.f. Section lHlB)l . 
and the spin-spin relaxation rates in 0.6% Gd doped EuO 
are longer than in pure EuO. 

The addition of 0.6% Gd in the EuO matrix gives rise 
to temperature dependent phenomena associated with a 
static magnetic inhomogeneity. The temperature depen- 
dent broadening mechanism initiates at about 30 K, the 
temperature at which the conductivity starts decreasing 
dramatically according to the results of Godart et alm^ 
and Samokhvalov et alH This suggests that the phe- 
nomena we observed at the level of the hyperfine field 
are likely to be related to the transport properties of the 
material and in particular to CMR. 

Magnetic entities such as magnetic polarons are 
thought to play a crucial role in these materials. In par- 
ticular, they are thought to be the main cause of CMR. 
We will show below that a simple model assuming that 
magnetic entities are present in the material and give rise 
to the static NMR line broadening, which we observe at 
higher temperatures, is not consistent with the observed 
relaxation times. 
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Suppose that in order to explain the narrow lines at 
low temperatures, we assume the magnetic entities dif- 
fuse rapidly at low temperatures, giving rise to motional 
narrowing. Thus they would create a fluctuating field at 
the nuclear site with correlation time r. The relaxation 
time deriving from the presence of this fluctuating field 
could then be written 

1 _ Slu 2 t 

¥2 ~ 1 + (Slut) 2 ( > 

where 8lo is the distribution in NMR frequency caused 
by the fluctuating magnetic fields Since we observed 
that the line is broadened at high temperature, we as- 
sume that at high temperature the correlation time is 
long, i.e. T2 = t. This is consistent with the model of 
magnetic polarons, since magnetic polarons are thought 
to be trapped at high temperature (leading to a decrease 
in conductivity) . At low temperature we assume that the 
entities are moving, which means that r is short. In this 
case, I/T2 = Slo 2 t and the line is motionaly narrowed. In 
between the high temperature regime and the low tem- 
perature regime, we expect that Slot = 1 at a given tem- 
perature. At that temperature, I/T2 is maximum and 
T2 = 2r = 2/Slu. Assuming that Slo is temperature in- 
dependent, that means that the relaxation rate is about 
equal to the static linewidth at high temperatures. But 
in fact, at all temperatures I/T2 is much less than the 
observed line width. Thus we cannot explain the data 
consistently with this model. Instead, we can argue from 
NMR that the material at low temperatures is spatially 
uniform, but as the temperature rises, non- uniformities 
set in. 

We also measured the effect of an external magnetic 
field of 4T on the temperature dependence of the FWHM 
of the line. A field of 4T was strong enough to saturate 
the magnetization, and as a consequence a non-negligible 
demagnetizing field was present in the non-perfectly el- 
lipsoidal sample. Because of the inhomogeneity in the 
demagnetizing field, we measured a FWHM of about 
0.5 MHz at 10 K, which was about 2.5 times broader than 
the zero-field FWHM at the same temperature. However, 
the FWHM did not increase with increasing temperature, 
and it was still about 0.5 MHz at 50 K. We concluded that 
the temperature dependent static inhomogeneities that 
we measured at zero-field were not present when there 
was a strong external magnetic field. This observation 
is an additional suggestion that the broadening might be 
linked to transport properties of the material. Indeed, 
since the conductivity of CMR materials is increased by 
the application of an external magnetic field, applying 
a field might perhaps be viewed as equivalent to lower- 
ing the temperature, as far as the transport properties 
are concerned^! Therefore, at a given temperature we 
can expect the line to be narrower in the presence of an 
external field than in zero-field. 

The center of the 2% and 4.3% Gd doped EuO 
lines are considerably shifted towards higher frequency 
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FIG. 9: Zero-field lineshape of Eu 153 in 2% Gd doped EuO vs. 
temperature. The amplitude of the lineshape was multiplied 
by the coefficient shown on the left side of each curve. 

(140.6MHz and 140.0MHz respectively). However, we 
did not observe a monotonic increase of the shift with 
the doping, the shift being smaller for 4.3% Gd doped 
EuO than for 2% Gd doped EuO. We also observed that 
increasing the Gd doping strongly influences the shape 
of the 153 Eu line (c.f. Fig.^J). In particular, the lines 
of 2% and 4.3% Gd doped samples were considerably 
broader than in the case of pure EuO and 0.6% Gd doped 
EuO and we could not distinguish the central transition 
from the other transitions. In order to decide whether 
or not this broadening mechanism was magnetic in ori- 
gin, we measured the 151 Eu lineshape of the 4.3% sample 
and we compare the 153 Eu and 151 Eu lineshapes. The 
frequency of the 151 Eu lineshape was multiplied by the 
ratio 7n/ 7n- Since the two curves had a similar 
shape, we concluded that the broadening mechanism was 
of magnetic origin. We do not know exactly what inter- 
action causes this large broadening, but we do know from 
the T2 measurements (c.f. Sect. lIIIBjl that the broaden- 
ing is a static broadening. In conclusion, we observed a 
remarkable difference between the lineshape of samples 
with none or low Gd doping and the lineshape of samples 
with higher Gd doping. 

We also measured the temperature dependence of the 
153 Eu lineshape in 2% and 4.3% Gd doped EuO. We plot- 
ted the measurements for 2% Gd doped EuO in Fig.|5] A 
similar temperature behavior was observed in 4.3% Gd 
doped EuO. We performed point- by-point measurements 
to determine the lineshapes. Note that at temperatures 
higher than about 80 K we could not detect the reso- 
nance. This failure was probably due to the fact that the 
line was too broad and therefore the signal was too weak. 

As shown in Fig. [51 the lineshape becomes broader with 
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FIG. 10: FWHM of the zero-field lineshape of Eu 153 in pure 
EuO and 0.6%, 2% and 4.3% Gd doped EuO vs. temperature. 



Gd doped) samples, and samples with higher Gd concen- 
tration both in static and dynamic magnetic properties. 
Below are the main outcomes of this study. 



1. Exchange integral J vs. doping 

We observed that the relaxation times in pure and 
Gd doped EuO were in good agreement with a law de- 
rived from spin- wave theory for temperatures as high as 
T = 0.6 Tc- From the temperature dependence of the 
relaxation times, we inferred the value of the exchange 
integral J as a function of Gd doping. We observed an 
abrupt increase of the amplitude of the exchange interac- 
tion with doping. This is consistent with the dependence 
on Gd concentration of the Curie temperature for which 
different theoretical explanations have been proposed by 
Mauger et aZ.jSLffi Nolting and Oles^tt and more recently 
Santos and Noltingi^ 



increasing temperature. Also, the structure of the line- 
shape becomes more and more complex with increasing 
temperature. In Fig.^H we plotted the FWHM of the 
lineshape as a function of temperature for 2% and 4.3% 
Gd doped EuO along with the temperature dependence 
of the FWHM of the lineshape of pure EuO and 0.6% 
Gd doped EuO. At temperatures above about 50 K the 
FWHM of 2% and 4.3% Gd doped EuO is ill-defined since 
the line has a complex shape. We determined that the 
temperature dependent broadening was magnetic in ori- 
gin and since 1 /T2 is by far smaller than the linewidth of 
the lines we measured, we concluded that the tempera- 
ture dependent broadening is a static magnetic broaden- 
ing. 

A possible explanation for the temperature dependent 
broadening in 2% and 4.3% Gd doped EuO is that the 
value of the exchange integral is distributed throughout 
the samples, i.e., different part of the samples have dif- 
ferent values of J. For any reasonable model of M(T), 
it would take a large change in Tc, hence of J, to ac- 
count for the spread in frequency we observed. However, 
a slight variation of J strongly influences the value of T\ 
because T\ oc J 9 / 2 . We observed that the measured Ti 
was, to a good approximation, independent of the posi- 
tion on the line. Therefore, we must conclude that there 
is a negligible distribution in J. 

V. CONCLUSIONS 

The study of pure and Gd doped EuO was motivated 
by the unconventional electrical and magnetic proper- 
ties observed in this system. Our NMR measurements 
led to results giving new information on the microscopic 
magnetic properties of Gd doped EuO. We observed a 
dramatic difference between pure or nearly pure (0.6% 



2. Variation of static magnetic broadening with 
temperature in Gd doped EuO 

We discovered that the lineshape of the Eu resonance 
in Gd doped EuO was broader than in pure EuO and 
that the broadening increases abruptly with Gd concen- 
tration. We also observed that the broadening increases 
with increasing temperature. We showed that the broad- 
ening was due to inhomogeneities in the local magnetic 
field acting on the Eu nuclei. We confirmed that the in- 
crease in broadening was not due to temperature depen- 
dent fluctuations of the electron spins. I/T2 was shown 
to be determined by these fluctuations, but was consid- 
erably smaller than the width of the broadening. Thus, 
the magnetic inhomogeneities were static at all temper- 
atures, at least at the time scale of NMR, which is of 
the order of several microseconds. Cooper et al., who 
studied Gd doped EuO by Raman spectroscopy, also ob- 
served magnetic inhomogeneities in Gd doped EuOiSii 6 ^ 
Our results confirm theirs but we can go further by af- 
firming that the inhomogeneities are static. 

In 0.6% Gd doped EuO, we observed that above about 
30 K the magnetic inhomogeneity increases rapidly with 
increasing temperature. According to previous trans- 
port measurements, the resistivity of samples containing 
a similar Gd concentration increases dramatically above 
about 30 This suggests that the broadening mech- 

anism is linked to the change in transport properties of 
the sample. Numerous models explain this change in 
transport properties by the formation of bound magnetic 
polarons. We could rule out the picture according to 
which the localized magnetic entities were highly mobile 
at low temperatures, giving motionally narrowed lines 
and high conductivity. 

In samples with higher Gd concentration (2% and 
4.3%), the observed static broadening is much larger than 
in 0.6% Gd doped EuO. Also, above about 30 K, some 
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structure develops in the lineshape along with the in- 
crease of magnetic inhomogeneity in the samples. We 
showed that this phenomenon was not due to a distribu- 
tion in the value of the exchange integral J. 

We must conclude that at low temperatures the mag- 
netic phase of the Gd doped samples is fairly homo- 
geneous, and that static magnetic inhomogeneities are 
formed when the temperature is increased. From our re- 
sults, we deduce that a theory to be valid has to include 
a spatially homogeneous and temperature independent 
J, but include the presence of magnetic inhomogeneities. 
Therefore, the explanation of the magnetization curve 
given by Borukhovich et al. is incompatible with our re- 
sults since it assumes a distribution of J in the sampled 
Although both the theories of Mauger et alS^ and Nolt- 
ing et al£SL assume a spatially homogeneous J and seem 
to correctly describe the doping dependence of the Curie 
temperature, they both suppose that J is temperature 
dependent, which in the light of our results cannot be the 
case. In addition, these models assume that the magne- 
tization is uniform within the sample, i.e. they do not 
take into account any magnetic inhomogeneities. 

There remains the possibility that the charge carrier 
trapping above about 30 K that accounts for the dra- 
matic temperature dependence of the electrical resistiv- 
ity causes a distribution in hyperfine field, presumably 
because the charge carriers are not all in the same spin 
state. This description is compatible with a theory based 



on the formation of magnetic polarons taken to be static 
on the NMR time scale. However, to the best of our 
knowledge, a complete model containing all the key fea- 
tures of the Gd doped EuO systems does not yet exists. 
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APPENDIX A: THREE-MAGNON RELAXATION 
PROCESS 

Below we present the derivation of the formula giv- 
ing the relaxation rate due to three-magnon relaxation 
process for the case of zero-field NMR in a ferromagnet 
with low anisotropy. Using a similar derivation than the 
one leading to (14) in the paper of Barak et a/.,— we 
obtained: 



1_ _ TrA 2 a 9 7 7 exp((e ki + gf i B \H\)/k B T) 

T x ~ 8HS J J exp((e fci + gf i B \H\)/k B T) - 1 
o o 1 

cx P(( e fc 2 +9VB\H\)/k B T) g(e ki )g{e ki )g(e k2 + ^k 2 ) de k 1 de k 2 
X exp((e k Jk B T + 9l i B \H\)) - 1 exp((e fci + e k ^ + 2 9f i B \H\)/k B T) - 1 ' 



where a is the lattice constant of the crystal (we used the 
fact that EuO and Gd doped EuO are cubic crystals), g 
is the g- factor, [i B is the Bohr magneton, H = H an + 
Ho + Hd m , where H an is the anisotropy field, Hq is the 
applied external field, and Hdm is the demagnetization 
field, e k is the spin- wave energy, and g{e k .) is the density 



of states of the spin waves. We then used the fact that 
the low energy spin waves give the principal part of the 
integral in (|Aip 4^ and we therefore determined g{e k ) by 

using the approximation e k = 2JSk 2 a 2 . We obtained: 



11 A 2 1 f f exp((e ki +g f i B \H\)/k B T) 



Ti 16(2tt)5 2 JS ■ HS (2 JS) 7 / 2 J J exp((e fr + gfi B \H\)/k B T) - 1 

o o 1 



eM(t k2 +WB\H\)/k B T) V 4x^+4^1 

X exp(( efc2 /fc B T + « ?MB |Jf|))-l exp(( efci + + 2g» B \H\)/k B T) -1 ^ ^ ( > 



r 



Beeman and Pincus used the fact that k B T ^> gfj, B \H\, 7.6 (k B T) 7 / 2 (c.f. (2.22) in the article of Beeman and 
and they found the integral in ljA2(l to be equal to 
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PincugSi)4ii We recomputed the integral in order to 
obtain a value adapted to the case of zero-field NMR 
measurements of Eu nuclei in EuO and Gd doped 
EuO. In that case, \H\ = \H an \ and at 4.2 K the ra- 
tio gnB\H an \/kBT is approximately equal to 8 • 10~ 3 . 
Since we observed that the three-magnon process is 
the main relaxation process for temperatures between 
about 15 K and 50 K, we had to take a substantially 
smaller value of g^B\H\/ksT . Taking T = 30 K and 
\H an (T = 30 K)\/\H an (T = 4.2K)| Si 2.5^ we ob- 
tained g[i B \H\jk B T = 4- 10~ 4 . With this value, we eval- 
uated the integral to be approximatively 10.5 (fcsT) 7 / 2 . 
This value is very close to the value obtained if we as- 



sume gfislHl/ksT = 0, and since there is no apparent 
reason to compute the integral with T = 30 K instead of 
any other temperature in the range 15 K< T <50K, we 
decided to use the result corresponding to the approxi- 
mation gfi B \H\/k B T = 0, that is 11.29 (k B T) 7 / 2 . This 
result leads to the expression given in 10. 

Since we were interested in evaluating J using l|A2J) . 
we had to determine the influence of a variation in the 
value of the integral on J. We calculated that if the 
factor 11.29 is replaced by the factor 10.5 obtained for 
T = 30 K, the value of J would be changed by about 1%. 
Therefore we decided that (JJJ could be used to determine 
a fairly accurate value of J from our data. 
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